Non-perturbative J p d model and ferromagnetism in dilute magnets 
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We calculate magnetic couplings in the J p d model for dilute magnets, in order both to identify 
the relevant parameters which control ferromagnetism and also to bridge the gap between first 
principle calculations and model approaches. The magnetic exchange interactions are calculated 
non-perturbatively and disorder in the configuration of impurities is treated exacly, allowing us 
to test the validity of effective medium theories. Results differ qualitatively from those of weak 
coupling. In contrast to mean field theory, increasing J p d may not favor high Curie temperatures: 
Tc scales primarily with the bandwidth. High temperature ferromagnetism at small dilutions is 
associated with resonant structure in the p-band. Comparison to diluted magnetic semiconductors 
indicate that Ga(Mn)As has such a resonant structure and thus this material is already close to 
optimality. 



is the coupling between d-orbital moments and the pri- 
marily p-orbital conduction carriers, x the concentration 
of d-orbitals and n the density of carriers. This estimate 
has even been influential in the search for new materi- 
als: attempts to increase T c at low doping have focussed 
on seeking large J p d and small bandwidth (large effec- 
tive mass). We shall show that the widely used Tq 
is in fact very misleading in the non-perturbative limit. 
Apart from potential applications Q , these materials are 
of fundamental interest for exploring ferromagnetism in 
a regime where geometric disorder is strong (because of 
the low concentrations of dopants) and the values of J p d 
are strong. Unfortunately these issues of principle have 
been obscured by insistence on the weak-coupling pic- 
ture, and has hidden richer underlying aspects. We argue 
that these systems are particularly interesting because 
the ferromagnetism, can only be understood if the crucial 
coupling J p d is treated non-perturbatively. An approach 
taking ab initio estimates of couplings and reliable treat- 
ment of the thermodynamics |5| can be successful, but 
the effective magnetic interactions found, which include 
a multitude of effects, have not been understood in sim- 
ple terms. Calculations must be taken for each individ- 
ual composition. A more systematic approach for finding 
promising new materials is possible if we have simple but 
accurate model systems. We can then identify which pa- 
rameters are important: J p d, the carrier density or the 
bandwidth of the conducting band? We can also examine 
assumptions of the ab initio approach: e.g. the Coherent 
Potential Approximation (CPA) for carrier disorder. The 
method we shall employ makes effective medium approxi- 
mations neither in the deriving of the magnetic couplings 
nor in calculating magnetic properties. Thus an aim is to 
make a bridge between ab initio and conceptually sim- 



The basis of our understanding of magnetism in 
metals is the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction Q. The magnetic exchange results from in- 
teraction of inner d- or f- localized orbitals and the outer 
s- or p- that form the conduction bands. In the standard 
form it is calculated with perturbative treatment of the 
coupling between the moments and the carriers on a def- 
inite Fermi surface, which is itself weakly affected by the 
presence of the magnetic moments. This leads to success- 
ful interpretation of oscillating exchanges in, for example, 
transition or rare earth metals. However the same con- 
cepts are often applied in situations where this perturba- 
tive approach is not really applicable. There are many 
situations where the coupling between moments and itin- 
erant carriers is not weak and it is important to determine 
changes in magnetic exchange in this regime, for example 
in the field of manganitesQ, double perovskites[3j or in 
materials modelled by Kondo lattices In another field 
of great current interest, diluted magnetic semiconduc- 
tors (DMS), a belief, incorrect in our view|5j, in the ap- 
plicability of weak coupling has resulted from an unhappy 
juxtaposition of two approximations: a simplified RKKY 
view of magnetic interactions and over-simplified molecu- 
lar field theory for the thermodynamics. This leads to the 
estimate for the Curie temperature T(f F — J^xn 1 ' 3 /t, 
where t — 1/m* is the inverse of the effective mass. J p d 
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pier model approach. A practical question is whether 
Tc can be increased by choice of other chemical dopants. 
In lightly doped III-V semiconductors, it will appear that 
Ga(Mn)As is in fact close to optimal. Here we discuss the 
theory appropriate to dilute magnets where the dopants 
provide well defined moments, but the concepts are rele- 
vant to a wider class of novel magnets ( magnets) in 
which there is no clear distinction between localized mag- 
netic moments and the itinerant carriers of the (doped) 
host HQ. 

The J p d model (or dilute Kondo lattice) is of a band 
of atomic d- orbitals coupled to an itinerant "p-band" of 
carriers via a Hund's rule coupling. In the manganites, 
for example, the "p" band is in fact be constructed from 
e g orbitals. 

H = - ^2 Ujc\ a c j(J + J& ■ s t + Vicl^Ci,, (1) 

tij = t for i and j nearest neighbors and zero otherwise. 
In the exchange between localized impurities spin and 
itinerant electron gas Ji is a random variable: Ji = J pc i 
and Vi = V if the site i is occupied by a magnetic impu- 
rity, and zero otherwise. Vi is here to mimic the effect of 
chemical substitution that accompanies the presence of 
the magnetic moment. Sj is the magnetic impurity spin 
operator at site i and Sj = cj Q (l/2<r Q( a)ci/3 the spin oper- 
ator, at the same site, of the itinerant electron gas. The 
d spins are distributed randomly at low concentration 
i on a regular lattice, for simplicity taken to be simple 
cubic. We remark that this model can be extended triv- 
ially to include several bands, as would be appropriate 
for models of double perovskites or manganites, where 
both the hopping and on-site potentials Vi may be 
random^J [ij. The results that will follow are non- 
perturbative in the J p d coupling and exact in the disor- 
der. A crucial effect that will be seen, and is neglected to 
low order in perturbation theory, is that as J p d increases, 
the conduction band is increasingly affected by the pres- 
ence of the magnetic ions. The calculations are made in 
the ferromagnetic phase, whose stability is tested a poste- 
riori, and therefore start from the configuration of spins 
S that are fully aligned in some direction S* = S for 
all occupied i. We can diagonalize in the Hilbert space 
(spins plus electrons) defined by the fully polarized d- 
spins provided that we ignore the transverse fluctuations 
of those spins, of order ^H3- Thus our theory is exact 
in the limit of large S — > oo For the electronic degrees of 
freedom this leaves a Hamiltonian quadratic in fermion 

operators. The effect of the d-orbitals is to produce an 

i J s 

effective spin-dependent term of Vf 1 — V±^j— on each 
occupied state, i.e. a random site potential for up (+)and 
down (-) carriers. For each configuration of disorder this 
Hamiltonian will be diagonalized exactly. 

In Fig^we show the single particle density of states for 
fixed JpdS = 5f and concentration x = 5% of magnetic 
impurities and different values of the d-orbital potential 




-6 -4 -2 2 4 6 8 -8 -6 -4 -2 2 4 6 8 10 




E/t 



FIG. 1: (Color online) Density of states. The bands of spin 
down and up are continuous and dashed respectively. Ener- 
gies are in units of t, J p dS = 5t and the concentration is 5%. 
The Fermi energy is indicated for three band fillings 7 

V/t = —1.2,-2.4,-3 and -5.4. The parameters were 
chosen to reproduce the general form seen in doped III- 
V semiconductors, as will be discussed in more detail 
below. For all calculations presented, the systems were 
sufficiently large (typically ~ 16 3 sites) that finite size 
effects are negligible. As V increases from the smallest 
value, a well defined impurity band for the minority band 
splits from the valence band. The Fermi energy is shown 
for different values of 7 = n^/x, where rih denotes the 
density of holes in a filled band. We will see below how 
the magnetic couplings are correlated with the form of 
the density of states and the position of the Fermi energy. 

The magnetic exchange interactions are derived at all 
distances between the different moments from the car- 
rier Green functions Gf,, fl^. This defines an effective 
(dilute) Heisenberg model for the magnetic moments, 
again valid for spins in the classical limit, H.Heis — 
— ■ JijSi ■ Sj. Exchanges are calculated by appl ying 
infinitesimal fields to the aligned ferromagnetic state[13j- 

" J —oo 

In this case Si = V^ — Vf~ . We note that the calculations 
use no effective medium approximation. In Fig. Gla) we 
show the couplings averaged over different impurity con- 
figurations. We vary the hole doping per d-orbital 7 for 
fixed x = 5% and V = —2 At, J p dS = 5t ( corresponding 
to FigOJb)). These clearly differ from RKKY for all fill- 
ings in that they do not oscillate around a mean value of 
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FIG. 2: (Color online) (a)Magnetic exchange as a function 
of distance, in units of lattice spacing, for different values 
of 7. Other parameters correspond to case (b) of Fig. (b) 
Exchanges for different values of the potential V at half filling. 
In inset, the exchanges rescaled: J(r)r 3 



zero and beyond 3 lattice spacings they essentially van- 
ish. While varying sharply with distance they do not 
change in sign, at least for small 7. The couplings at 
short distances increase monotonically with 7 but for suf- 
ficiently large 7 antiferromagnetic values appear at larger 
distances (see inset). In Fig. E{b) the J(r) with differ- 
ent potentials V are shown for fixed 7. As |V| increases 
the interactions increase at short distances but become of 
very short range. For \V\ small the longer range oscilla- 
tions appear clearly, as seen in the inset. Thus from FigEl 
there is a range of V and 7 where the couplings are all fer- 
romagnetic ( J(r) > 0). This is associated (see Fig. HO 3 )) 
with incipient development of a visible impurity band just 
at the band edge. This ferromagnetic "bias" results from 
the resonant form of the p-band as advanced in Ref . • 
In that paper we compared to exchan ge c alculated in a 
toy model between two resonant statesjl5j . arguing that 
particle-particle channels responsible for super-exchange 
should be suppressed. In contrast, in the present calcula- 
tion there is no artificial separation between the resonant 
states and the rest of the conduction band nor between 
different contributions (particle-hole and particle-particle 
channels). The ideas are coherent in that the J p d acting 
between the d- spin and the resonant p-state suppress 
particle-particle virtual states. For the largest values of 
\V\ the clear suppression of exchange for large distances 
(see inset to Fig. ETb)) is a precursor to pure magnetic 
double exchange type [la] ■ For J p d and V strictly infinite 
we expect the coupling to reduce to nearest neighbours: 
Jij — tij (cfcj)[l(| For JpdS finite, however, the exch ang e 
interactions extend further than nearest-neighbours |l7|. 
We now calculate the Curie temperatures Tc by the 



FIG. 3: (Color online) Curie temperature in units of t as a 
function of carrier density. In insert the same calculation for 
the RKKY couplings see ref Q Fig. 4 

self-consistent local random phase approximation (SC- 
LRPA)0 which treats geometric disorder exactly, as is 
particularly important in the dilute limit. SC-LRPA has 
been successful in reproducing experimental Tn of DMS 
and agrees with Monte Carlo simulations [l^, fl9| . Figs. El 
and 2] shows the Curie temperatures (Tc) as functions 
of 7 and V at half-filling respectively, and various val- 
ues of Jpd- Unlike the predictions of T^ F , T c increases 
monotonically neither with J p d nor with 7. For 7 = 0.1, 
for instance, Tc decreases with J p d\ Large values of J p d 
are primarily useful in broadening the range of stability 
of ferromagnetism with 7. For comparison we show in 
inset previous results in the RKKY regime 0|- At the 
smallest values of J p d the region of stability (note the 
difference in scale 7) approaches the RKKY results, as 
expected. In Fig. 0] we see that Tc is less sensitive to J p d 
in the large V limit. This is essentially the limit of dou- 
ble exchange. The ferromagnetic bias is essentially the 
physics of double exchange 20J, and the scale of Tc is t. 
Note we make the assumption that the ferromagnetic Tc 
calculated with couplings averaged over configurations is 
equal to that obtained, in the thermodynamic limit, with 
unaveraged couplings. We have checked this assumption 
and it is accurate. 

We can use our results as a simple model for doped 
Ga(Mn)As. We take parameters from photoemission |2 lj 
J pd « l.2eV,S = 5/2, J pd S « ieV. From the band- 
structure of GaAs: m* ps 0.5m e , lattice spacing oq = 
5.65A, this gives t w 0.7eV and thus J pd = At. The 
value of V/t = —2AeV was chosen to mimic the calcu- 
lated density of states of Ga(Mn)As. For these parame- 
ters and x = 5%, from Figure. E3 we obtain a maximum 
rpiband ^ at 7 = 0.3. Because there are two degen- 
erate hole bands and we have neglected all correlations 
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FIG. 4: (Color online) Curie temperature in units of t at half- 
filling as a function of potential V and varying coupling J p d- 
Points (a)-(d) correspond to the density of states of Figures 
l(a)-(d) respectively 



principle calculations for calculating magnetic exchanges, 
an approximation that has often been questioned. While 
we have presented results in the dilute regime, we note 
that our approach, which gives the exchange interactions 
in real space, can be used in more general situations at 
higher concentrations or with inhomogeneities, as may be 
applicable in other materials described by the J p d model 
and its extensions. For dilute magnets a region of param- 
eter space gives couplings that are all ferromagnetic. The 
reason^) is made explicit: optimality is associated with 
a resonant structure seen in the p-band. The potential 
V must be sufficiently strong that this occurs, but if it 
is too strong, couplings are too short-range to maintain 
long range order. Unlike previous mean field theories 
(Tq F ) : the temperature scale for ferromagnetism is con- 
trolled not by Jpd but by the bandwidth t. The physics 
is essentially that of double exchange in this regime. The 
values of Fermi energy, J p d and the potential shift V are 
important in obtaining optimal couplings. 



in the p-bands, we can simply extend to 2 independent 
bands: T*> and = 2T^ and « 130K, close to the measured 
value in annealed samples of Ga.9s(Mn)o.o5As. This is a 
crude comparison as the lattice is face-centred cubic and 
not cubic and, at first sight, the natural concentration to 
chose for uncompensated samples might be two half-filled 
bands 7 = 0.5 to give one carrier per dopant. We argue, 
however that 7 = 0.3 is comparable to uncompensated 
Ga(Mn)As since it is at this doping that the density of 
states (calculated from band structure j^) most resem- 
bles that of Figur^J), with Ep slightly above the maxi- 
mum in the impurity peak. The apparently remarkable 
agreement for such a simple model is found for the range 
of concentrations relevant to experiments. We conclude 
if Jpd is treated non-perturbatively we do not need exten- 
sions such as the 6-band Luttinger model,?]] to give ac- 
curate quantitative comparison. We can also understood 
from the position of the impurity band why Ga(Mn)As 
has a larger Tc than Ga(Mn)N and In(Mn)As which re- 
semble more cases (d) and (a) of Fig. ^ respectively 0. 
Since the estimated J p d and bandwidths for these mate- 
rials are all very similar, the important difference is the 
potential V. From the calculated Tqs shown in the Fig- 
ures and comparison with the density of states in Fig^ 
(b),we see that Ga(Mn)As has close to optimal parame- 
ters. This may explain failure to find higher values of Tq 
in lll-V semiconductors without invoking other effects: 
local clustering or increasing bandwidth. 

To conclude, we have shown that in the non- 
perturbative regime the J p d model has effective magnetic 
couplings very different from those of RKKY and resem- 
ble those of ab initio calculations for the DMS. This sim- 
ilarity supports the accuracy of CPA as used in the first 
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